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A simple new binary test for chaos has been proposed by Gottwald and Melbourne. We apply 
this test successfully to the Henon-Heiles and Lorenz systems, demonstrating its applicability to 
conservative systems, as well as dissipative systems. The binary test is effective for highly chaotic 
Hamiltonian systems and orbits on a strange attractor and is particularly useful as a marker of the 
transition from regularity to chaos. However, we find it is not able to detect more subtle instances 
of transient chaos. 



PACS 05.45.-a,82.40.Bj 

I. INTRODUCTION 

A promising new binary test for chaos has been pro- 
posed by Gottwald and Melbourne [1]. The test returns 
for a non-chaotic system and 1 for a chaotic system. It 
does not offer any quantitative information and so can- 
not quantify the degree of irregularity in the way a Lya- 
punov exponent might. However, it can provide a quick 
and useful diagnostic for chaos. 

Consider a base dynamical system x = f(x) of any 
dimension with orbits x(i). The proposed method to 
determine if it is chaotic is based on a Euclidean E(2) 
group extension of the underlying base dynamics [3,4]. 
The dynamics is explicitly extended to include two new 
variables (p, q) defined through 



p = 0(x) cos(w i), 
q = 0(x) sin(wot), 



(1.1) 



where </>(x) is any observable of the base dynamics and 
luq ^ is an arbitrary constant frequency that is needed 
to damp off any linear growth that may be common to 
both non-chaotic and chaotic orbits. If the observable cj) 
is drawn from a non-chaotic system, then the (p, q) sub- 
space will be bounded. But if the observable 4> is drawn 
from a chaotic system, then the motion in the (p, q) sub- 
space will be Brownian and unbounded. The extended 
variables diffuse through the subspace as the observables 
of the base dynamics jump around unpredictably. 

Gottwald and Melbourne introduced a binary test for 
chaos by defining a mean-square displacement 

M(t)= lim I f (p(t + r)-p(T)) 2 dr (1.2) 

and characterizing the behaviour of M(t) through 

log M(t) 



if = lim 



{-too logi 



(1.3) 



If there is no chaos then the motion in p{t) is bounded 
and if — > 0. If the base dynamics is chaotic then p(t) will 
exhibit Brownian diffusion so that Ap — > t 1 ^ 2 as t — > oo 
and K — * 1. The if -test can be used for both continuous 
and discrete dynamical systems even when the precise 
underlying base dynamics is unknown [1]. 



This test has certain advantages of simplicity and was 
shown to confirm the results of testing for chaos using 
Lyapunov exponents for the forced van-der-Pol system 
[1] . We show here that the if -test also confirms the trend 
of the Lyapunov exponents for the simple Hamiltonian 
system of Henon and Heiles [2] . So the if -test is appli- 
cable to Hamiltonian as well as dissipative systems. 

While the test works well as a diagnostic of the tran- 
sition from regularity to chaos we express some reserva- 
tions about the use of the test on an orbit by orbit basis. 
Most importantly, we found the method difficult to in- 
terpret for systems that experienced chaotic transients in 
contrast to the easy interpretation of dissipative systems 
near a strange attractor, or in highly chaotic Hamilto- 
nian systems. To illustrate this we use the Lorenz model 
[5] to show that transient chaos, which can be found in 
the Lorenz system by other means, goes undetected by 
the if -test. Orbits that enter highly chaotic regions of 
phase space and then depart into regular regions are not 
Anosov and may not show the diffusive motion in phase 
space required to yield a if — > 1 result from the test. 
These transient systems may not satisfy the criteria for a 
system within the remit of the if -test. This signals one 
of the possible limitations of the K-test. 



II. HENON-HEILES 

The Henon-Heiles system is an ideal Hamiltonian sys- 
tem on which to evaluate the utility of the if -test. It pro- 
vides a simple and well understood model for the motion 
of stars in a galactic potential as well as for the motion 
of non- linearly-coupled molecules. The transition from 
regularity to chaos as the energy is increased is well doc- 
umented and can be identified with Poincare surfaces of 
section, Lyapunov exponents, and the Painleve property. 

The Henon-Heiles Hamiltonian is 

H=^(pl+p 2 y + x 2 +y 2 )+x 2 y~±y 3 (2.1) 
with equations of motion 

X = p x 

V=Py 
p x = -x - 2xy 
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p y = -y-x 2 +y 2 . (2.2) 

We add the Euclidean extension 

p = xcos(6) 
q = x sin(#) 

9 = uj (2.3) 

where we have chosen <j> — x. We evolve these equations 
numerically and check that the Hamiltonian remains con- 
served throughout the simulation. Our results are sum- 
marized in Figs. 1 and 2. The orbits selected for display 
have identical initial conditions up to the energy H which 
was varied from orbit to orbit. The value of K(t) is mea- 
sured as the slope of the line ln(M(t)) versus ln(t) and 
is plotted in Fig. 1 as a function of the energy H. The 
transition from non-chaotic orbits (K — > 0) to chaotic 
orbits (K — > 1) as the energy increases is clearly demon- 
strated. The transition confirms the break up of tori 
in the phase space for the base dynamics (and therefore 
the non-integrability of the base dynamics) as the energy 
grows. We also confirm the onset of chaotic motion for 
these orbits with a numerical determination of the prin- 
cipal Lyapunov exponent shown as a function of H in 
Fig. 2. 



A few comments on Fig. 1 should be made. The con- 
vergence of the value of K towards either or 1 can 
be improved if the numerical simulations are allowed to 
run far longer and a more precise extrapolation method 
is used so as to avoid any early transients. Of course, 
this renders the application more numerically intensive 
and slower without yielding a very great improvement in 
convergence. There is also a clear dependence of the mo- 
tion in the Euclidean subspace on the specific value of uio 
employed. However, the overall trend in K is typically 
the same for various values of uj - 

While Fig. 1 does show the required transition from 
to 1, it is important to stress how the conservative sys- 
tem differs from the dissipative system. For a dissipative 
system all trajectories, regardless of initial conditions, 
are drawn onto the same attractor. Therefore the K-test 
should show no dependence on initial conditions. This 
is not true for a Hamiltonian system. There can be a 
mixture of regular and irregular orbits and so the test 
result can depend on initial conditions. Only for a com- 
pletely ergodic system will K — > 1 independent of initial 
conditions. Also, since chaotic transients can often arise 
in conservative systems, the K-test may not always give 
such a crisp transition. Chaotic transients in dissipative 
systems are discussed further in the next section. 




FIG. 1. The asymptotic value of K versus the energy H 
for the Henon-Heiles system. 




III. LORENZ 

The standard Lorenz system demonstrates a transition 
from no chaos to transient chaos and then further to a 
full chaotic attractor. The standard system is 

X = -a(X-Y), 

Y = -XZ + rX -Y, 

Z = XY-bZ, (3.1) 

with a = 10 and b = 8/3 and we vary the constant r. 
The transition to chaos is summarized in Ref. [6]. For 
r < 1 there is one fixed point &tX = Y = Z = and 
no chaos. For 1 < r < 13.96 there are two additional 
attractors, one at 

X R = (b(r-l)) 1 ' 2 
T^(&(r-l)) 1/2 
Z R = r-l, 



(3.2) 



FIG. 2. The asymptotic value of the principal Lyapunov 
exponent versus the energy H for the Henon-Heiles system. 



and another at (X L ,Y L , Z L ) = (—X R , —Y R , Z R ), but still 
no chaos. A non-chaotic orbit drawn onto the simple 
fixed point (X R , Y R , Z R ) with (X, Y, Z) — > is shown in 
Fig. 3. A closely neighbouring orbit will follow a similar 
path onto the same attractor since there is no exponen- 
tial sensitivity to initial conditions, no positive Lyapunov 
exponent, and no chaos. 

However, for values of r > 13.96 a non-attracting form 
of transient chaos develops around a homoclinic orbit in 
phase space, first noted by Kaplan and Yorke [7]. Two 
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orbits undergoing transient chaos are shown in Fig. 4. 
One orbit winds around a bit before being drawn onto 
the regular attracting point (Xl,Yl, Zi) while its very 
near neighbour ends up on (Xr, Yr, Zr). 





FIG. 3. A nonchaotic orbit for r = 10 being drawn onto 
the regular attracting fixed point at (Xr, Yr, Zr). The orbit 
began with initial conditions X(0) = 1.99895, Y(0) = 0, and 
Z(Q) = 9.99005. 
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FIG. 4. An orbit in the (X, Z) plane with r = 18 exhibit- 
ing transient chaos before gliding towards the (not strange) 
attractor at (Xl, Yl, Zl)- Left: The orbit began with initial 
conditions X{0) = 1.99885, Y{0) = 0, and Z(0) = 9.99005. 
Right: A neighbouring orbit with X(0) = 1.99895, Y(0) = 0, 
and Z(0) = 9.99005 is drawn onto the other regular attractor 
at (Xr,Yr,Zr). 

We choose to use the method of fractal basin bound- 
aries to locate regions of transient chaos. An initial slice 
in phase space is color-coded according to whether the 
orbit eventually lands on the fixed point (Xr, Yr, Zr) 
(black) or the fixed point (X L ,Y L , Z L ) (white). In the ab- 
sence of chaos all the basins will have smooth boundaries. 
In the presence of chaos the boundaries become fractal, 
demonstrating both extreme sensitivity to initial condi- 
tions and a chaotic mixing of orbits. Such fractal basin 
boundaries are illustrated in Fig. 5; the orbits shown in 
Fig. 4 are drawn from along the fractal. They have nearly 
identical initial values but divergent outcomes, a charac- 
teristic of chaotic transients. 




FIG. 5. Upper Figure: A fractal basin boundary in (X, Z) 
for r = 18. The range of initial conditions is —10 < X(0) < 10 
and < Z(Q) < 10 while for all orbits Y(0) = 0. Middle 
Figure: A detail of the fractal basin boundary in (X,Z). 
The range of initial conditions is 1 < X(0) < 2 and 
9.9 < Z(Q) < 10 while again Y(Q) = 0. Lower Figure: 
A detail of the detail. The range of initial conditions is 
1.99 < X(0) < 2 and 9.99 < Z(0) < 10 with Y(0) = 0. 




FIG. 6. An orbit in the (X, Z) plane along the strange 
attractor for r — 30. The orbit began with initial conditions 
X(Q) = 1.99895, Y(0) = 0, and Z(0) = 9.99005. 

For r > 24.06 the transient chaos gives way to chaos on 
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attractors which merge into the famous Lorenz strange 
attractor beyond r > 24.74 [6]. An orbit that drifts onto 
the strange Lorenz attractor is shown in Fig. 6. 

We apply the K-test to the Lorenz system and show 
that it effectively marks the transition from non-chaotic 
motion to chaos on a strange attractor at r ~ 24.74 where 
we see K rising from to 1. However, the test is unable 
to pick up the chaotic transient behaviour for values of 
13.96 > r > 24.06. Notice that a similar random scan 
through the Lyapunov exponents in Fig. 8 also misses the 
transient episodes although short time exponents along 
the fractal basin boundary can be isolated. 



subtle chaotic transient in the early motion. 
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FIG. 7. The asymptotic value of K versus the parameter r 
for orbits with inital values X(Q) = 1.99895, Y(0) = 0, and 
Z(0) = 9.99005. 




FIG. 8. The principal Lyapunov exponent versus the pa- 
rameter r for orbits with inital values X(0) = 1.99895, 
Y(0) = 0, and Z(0) = 9.99005. 

Specific orbits of Fig. 7 can be isolated to illustrate the 
behaviour in the (p, q) subspace explicitly. The orbit at 
r = 18 corresponds to the right-most orbit in Fig. 4. The 
orbit winds around before drifting onto the fixed point 
(Xr, Yr, Zr). This is reflected in the (p, q) subspace by 
the stray steps taken before the orbit moves onto the 
smooth, bounded ring of the regular attractor as shown 
in Fig. 9. This results in a regular oscillation in the mean- 
square displacement M(t) that will eventually die away 
to give K — > 0, as shown in Fig. 7. The if -test reflects 
the regularity of the attractor and is unable to detect the 




FIG. 9. The bounded regularity of the (p, g)-extension for 
an orbit sampled from Fig. 7. The orbit corresponds to 
r = 18 with inital values X(0) = 1.99895, Y(0) = 0, and 
Z(0) = 9.99005 (see Fig. 4). 

By contrast, an orbit at r = 30 in Fig. 7 is drawn 
onto the Lorenz strange attractor and its strongly chaotic 
behaviour is detected by the if -test. The chaotic motion 
is well reflected by the Brownian diffusion evidenced in 
the (p, q) subspace shown in Fig. 10. This orbit leads to 
K — > 1 as expected. 




FIG. 10. The figure shows the unbounded Brownian like 
diffusion of the (p, g)-extension for a chaotic orbit for which 
K ~ 1 in Fig. 7. The orbit corresponds to r — 30 with inital 
values X(0) = 1.99895, Y(0) = 0, and Z(Q) = 9.99005 and 
moves onto the Lorenz strange attractor. 



IV. SUMMARY 

We have confirmed that the if -test provides a simple 
and easy diagnostic for the transition from regularity to 
chaos for Hamiltonian as well as dissipative systems. For 
the test to be effective, an orbit must spend sufficient 
time on the hyperbolic chaotic attractor in a dissipa- 
tive system, or be confined to a highly chaotic region 
of phase space if the system is Hamiltonian. However, 
for chaotic transients which move into a chaotic region of 
phase space and then out into a regular region of phase 
space the motion is not consistently Brownian. Conse- 
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qucntly, the test can yield ambiguous and confusing re- 
sults although a qualitative look over the (p, q) subspace 
can provide guidance as to the transient irregularity of 
the base dynamics. Additionally the K-test is not ideal 
for relativistic settings since it depends on the time coor- 
dinate used and so, like the Lyapunov exponents, is not a 
covariant indicator of chaos. These limitations of the K- 
test are not unexpected: no one probe of chaos can suit 
every scenario. Nor are they fatal to its utility. Rather, 
we hope that they will help to map out the territory over 
which this simple test is a reliable guide to the presence 
of chaos 
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